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ABSTRACT 

We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), 
HS1700.850.1, the brightest SMG (S H80;/tll =19.1 mJy) detected in the JCMT/SCUBA-2 Bary- 
onic Structure Survey, based on the detection of its 12 CO line emission. Using the IRAM 
PdBI-WIDEX with 3.6 GHz band width, we serendipitously detect an emission line at 
150.6 GHz. From a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, 
we confirm the identification of this line as 12 CO(5-4) at z = 2.816, meaning that it does 
not reside in the z ~ 2.30 proto-cluster in this held. Measurement of the 870/jm source 
size (<0.85") from the Sub-Millimeter Array (SMA) confirms a compact emission in a 
Ss 70 /rm =14.5 mJy, L/ /? ~ 10 13 L 0 component, suggesting an Eddington-limited starburst. 
We use the double-peaked 12 CO line profile measurements along with the SMA size con¬ 
straints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses 
of HS 1700.850.1. While HS 1700.850.1 is one of the most extreme galaxies known in the 
Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution 
in this held. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), 
and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively 
common, with the brightest SMGs rarely occupying the most overdense regions at z =2-4. 

Key words: galaxies: abundances - galaxies: high-redshift - submillimeter: galaxies. 


1 INTRODUCTION 

Blank-field mm and submm continuum surveys have discovered 
hundreds of dusty, star-forming submm galaxies (SMGs) over the 
past decade (e.g., Smail et al. 2002; Borys et al. 2003; Greve et al. 
2004; Coppin et al. 2006; Bertoldi et al. 2007; Scott et al. 2008; 
Weiss et al. 2009), while Herschel-SPIRE has mapped close to 10 3 
deg 2 in the 500, 350, 250/.;m bands, leading to ~10 4 SMGs be¬ 
ing detected (e.g., Oliver et al. 2010). At higher flux limits SPT 
has mapped 2500 deg 2 (e.g., Mocanu et al. 2013), and Planck has 
mapped the entire sky (Planck collaboration et al. 2014), finding 
many extreme and/or gravitationally lensed SMGs. The search for 
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the most luminous star forming galaxies from 100’s of deg 2 in 
the Universe has ensued (e.g,. Fu et al. 2013, Ivison et al. 2013). 
However some of the SMGs discovered in early ~100 arcmin 2 sur¬ 
veys remain the intrinsically brightest, most extreme known (e.g., 
SMM J02399 - Ivison et al. 1998, 2010, GN20 - Pope et al. 2006; 
Younger et al. 2008, and COSMOS-AzTECl - Scott et al. 2008; 
Younger et al. 2010). 

Determining the redshifts of unlensed SMGs has largely been 
a process of using weak counterparts in the rest-frame ultraviolet 
and optical to make spectroscopic determinations (e.g,. Chapman 
et al. 2003, 2005), or combining near- and mid-IR photometry to 
obtain reliable photometric redshift estimates (e.g., Wardlow et al. 
2011; Simpson et al. 2014). However, because the high extinction 
of SMGs means that they often have only very faint (if any) coun¬ 
terparts in the optical/near-IR bands, the spectroscopic and even 
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Figure 1. Left to right: Five cutout images 25" on a side of the longer wavelength data for HS 1700.850.1 (850/.im, 450//m, 24//m. 8//m. 4.5//m). highlighting 
the likely identification of the SMG with an IR-luminous component (black circle shows the 15" SCUBA-2 beam size). Second from right: A 9" x9" zoom 
in from the SMA compact+extended configuration 870//m map, with a ~ 0.74" synthesized beam with PdBI 12 CO(5-4) contours overlaid. Most of the 
SCUBA-2 flux (14.5 mjy) is recovered in a single compact component, unresolved with a 1" beamsize. No other significant sources are detected within the 
15" diameter of the SCUBA-2 beam. Further, the compact configuration data on its own measures a flux of 19±3 mjy, comparable to the SCUBA-2 flux, and 
suggestive that the source may be marginally resolved on scales > 1". Far right: 9" cutout of the HST -ACS F814W image centred on the SMA detection, 
with IRAC ch-2 contours overlaid, shows F1S1700.850.1 to be well identified with a faint tad-pole like I >\ /-> ~ 26 galaxy, which is marginally detected in 
ground based imaging (' R. = 26.1). 


photometric redshift distributions have remained incomplete and 
biased (e.g., Simpson et al. 2014). 

The largest spectroscopic SMG redshift survey to date was 
based on radio-identified SMGs (Chapman et al. 2005). The radio 
identification is sometimes inaccurate (Hodge et al. 2013) and may 
bias the redshift distribution since radio emission may remain unde¬ 
tected even in the deepest radio maps for sources at 2 > 3. Optical 
spectroscopic followup to ALMA-identified SMGs (Simpson et al. 
2014; A. Danielson et al. in prep) have removed the identification 
bias, yet cannot remove the spectroscopic bias. 

An alternative route to determine the redshift of an SMG is 
through observations of 12 CO emission lines at cm or mm wave¬ 
lengths. The 12 CO lines arise from the molecular gas, the fuel 
for star formation, and can thus be related unambiguously to the 
submm continuum source. Therefore, these observations do not re¬ 
quire any additional multi-wavelength identification and circum¬ 
vent many of the problems inherent to optical spectroscopy of 
SMGs. However, to date the sensitivity of mm-wave facilities has 
limited the approach to followup of bright, gravitationally lensed 
SMGs (e.g., Weiss et al. 2009; Swinbank et al. 2010; Harris et al. 
2012; Vieira et al. 2013; Zavala et al. 2015). 

While the IRAM Plateau de Bure Interferometer (PdBI) is 
sensitive enough to detect 12 CO in unlensed SMGs (e.g. Both- 
well et al. 2013), its 3.6 GHz bandwidth receivers limit this ap¬ 
proach due to the large time requirement for blind searches of 12 CO 
lines in redshift space via multiple frequency tunings. Walter et 
al. (2012) observed the HDF-N field with the PdBI over the full 
3mm window to a similar depth as the Bothwell et al. (2013) sur¬ 
vey, requiring 110 hrs on-source. In addition to blindly detecting 
two faint 12 CO emitting galaxies, they identified the redshift of the 
Ss 50 /im ~5mJy HDF850.1 as 2 = 5.3. The 8-GHz instantaneous, 
dual-polarization bandwidth of ALMA, coupled with its sensitiv¬ 
ity, improves this approach, but still makes it a relatively expen¬ 
sive prospect to blindly obtain 12 CO redshifts for even the brightest 
SMGs with Ssso/jm ~8mJy. 

In this paper, we describe a photometric redshift guided, mm- 
wave search for 12 CO line emission from HS1700.850.1, one of 
the brightest unlensed SMGs ever discovered, with S’ssop.m = 
19.1 mjy. The target field was chosen from the 24 survey fields 
of the Keck Baryonic Structure Survey (KBSS - e.g., Rudie et 
al. 2012). Several of these fields show strong galaxy overdensi¬ 
ties (e.g., Steidel et al. 2000, 2005, 2010), and together represent 
a complete, unbiased census of overdensities in a well-calibrated 
field galaxy survey at 2 > 2. The HS 1700+64 field in which 
HS1700.850.1 was found hosts a well characterized 2 = 2.30 pro¬ 


tocluster, studied in detail in several works (Steidel et al. 2005; 
Shapley et al. 2005; Digby-North et al. 2010). The faintness of 
HS1700.850.1 at near-IR/optical wavelengths precluded spectro¬ 
scopic observations, despite lying in this well studied field. 

We use cosmological parameters = 0.3, A = 0.7, and H 
= 70 km s- 1 Mpc -1 throughout the paper; at 2 = 2.82, this cor¬ 
responds to an angular scale of 7.94 kpc arcsec -1 , or 0.48 Mpc 
arcmin~ 1 . 


2 OBSERVATIONS AND STRATEGY 

The HS 1700+64 field was targeted with SCUBA-2 at the JCMT 
as part of a submm-wave followup of six KBSS fields (K. Lacaille 
et al. in prep), covering ~ 0.19 deg 2 . The brightest source in this 
850/rm survey is HS 1700.850.1 with S850//m=19.1±0.8rnJy. Fig. 
1 shows cutouts of HS 1700.850.1 at SCUBA-2 and Spitzer wave¬ 
lengths, revealing a likely identification at 4.5pm. 

The HS 1700+64 field contains a massive proto-cluster at 2 = 
2.30 (Steidel et al. 2005), and given its 1.5' proximity to the clus¬ 
ter centre, HS 1700.850.1 was considered a likely cluster mem¬ 
ber. HS 1700.850.1 has two candidate Spitzer -IRAC identifications 
within 4" of the 850^tm centroid, one of which is associated with a 
well detected complex of 2 = 2.305 galaxies. IRAM-PdBI obser¬ 
vations targeting 12 CO(3-2) and [Cl] at 2 = 2.305 in the 3mm and 
2mm bands respectively showed no detectable line emission, indi¬ 
cating this redshift was incorrect (Fig. 2). However, the WIDEX 
receiver’s 3.6 GHz wide spectrum revealed a strong line detected 
serendipitously at 150.6 GHz near the edge of the chosen 2mm- 
band setting (Figs. 2,3). Two followup observations targeted the 
three most likely redshifts, assuming this 150.6 GHz line was a 
12 CO rotational transition, as described in § 3.1. We also draw on 
Spitzer- IRAC, ground-based and HST optical datasets described by 
Steidel et al. (2011). All flux measurements are shown in Fig. 4, 
and discussed in § 3.1. 

2.1 JCMT/SCUBA-2 Observations 

SCUBA-2 observations were conducted in Band 2 weather con¬ 
ditions (T 225 GHZ < 0.08) over 6 nights between 23rd May and 
15th Oct. 2013, totalling 15 hours of on-sky integration in indi¬ 
vidual 30 min scans. A standard 3' diameter “daisy” mapping pat¬ 
tern was used, which keeps the pointing centre on one of the four 
SCUBA-2 sub-arrays at all times during exposure. Data reduction 
followed standard recipes (e.g., Geach et al. 2013), with individual 
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Figure 2. A composite of all spectral scans towards HS 1700.850.1, shown at a velocity resolution of 60 kms -1 (arbitrary flux scaling). 12 CO emission lines 
are seen at ~ 90.2 and 150.6 GHz, the latter corresponding to 12 CO(5—4). The originally targeted z = 2.31 12 CO and [Cl] line frequencies are highlighted 
(magenta), with a curious excess seen at the [Cl] frequency. The four tuning bands are marked with heavy lines, showing a slight overlap near 101 GHz. 
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Figure 3. Spectra of the 12 CO(3-2) and 12 CO(5-4) lines towards 
HS 1700.850.1. The spectral resolution is 60 km s -1 for both lines. We find 
only marginal evidence for any spatial separation in the two peaks of the 
CO lines < 1". Double Gaussian fits to each line profile are overlaid, em¬ 
phasizing the ~50% higher peak flux in the blue component of the CO(3-2) 
profile. Line ratio and dynamical modelling issues are discussed in detail in 
the text. 


Table 2. Photometry of HS 1700.850.1 


Band 

Flux 

Instrument 

91.7GHz 

0.25±0.07 (mJy) 

PdBI 

100.8GHz 

0.31 ±0.13 (mJy) 

PdBI 

103.8GHz 

0.38±0.08 (mJy) 

PdBI 

149.2GHz 

1.39±0.15 (mJy) 

PdBI 

850/im 

19.1 ±0.8 (mJy) 

SCUBA-2 

870/i,m 

14.5± 1.1 (mJy) 

SMA 

450/i,m 

45±6 (mJy) 

SCUBA-2 

24 pm 

171 ±7 (pjy) 

MIPS 

S.Ofim 

13.2± 1.2 (pjy) 

IRAC 

4.5 fim IRAC 

10.2±0.9 (pjy) 

IRAC 

-ff 3 (2.18pm) 

23.20±0.13 (AB mag) 

P200/WIRC 

J(1.25pm) 

24.4±0.2 (AB mag) 

P200/WIRC 

/(F814W) 

25.8±0.2 (AB mag) 

HST/ACS 

n a 

26.1 ±0.3 (AB mag) 

WHT/PFC 

9' 

27.0±0.3 (AB mag) 

WHT/PFC 

u' 

>27.1, 3cr (AB mag) 

WHT/PFC 


“ The 77. filter is centered at 6900A and thus redder than the standard R - 
close to mid-way between SDSS r' and i' in effective wavelength 


Table 1. PdBI observations log. 


Run-ID 

center freq 
(GHz) 

t? , 

int 

(hr) 

Depth 6 

(mJy) 

obs date 

beamsize 

S14ck001 

103.8 

3.34 

0.60 

09/06/14 

5.1" x4.2" 

S14ck002 

149.2 

3.45 

1.12 

11/06/14 

3.5" x2.7" 

D14d001 

100.8 

1.50 

0.95 

28/07/14 

O 

X 

’7}- 

D14d002 

91.7 

3.15 

0.43 

04/08/14 

4.8" x 4.1" 


a 5-Antenna integration. 
b per 50MHz channel. 


nal from each bolometer’s time stream is then re-gridded onto a 
map, with the contribution to a given pixel weighted according to 
its time-domain variance. Filtering of the time-series is performed 
in the frequency domain, with band-pass filters equivalent to angu¬ 
lar scales of 2 "< 6 < 120", along with IOct deviation spike re¬ 
moval in a moving boxcar, and DC step corrections. Finally, since 
we are interested in (generally faint) extragalactic point sources, we 
apply a beam matched filter to improve point source detectability, 
resulting in a map that is convolved with an estimate of the 850pm 
beam. The map reaches central noise levels of ers5o=0.51 mJy and 
0-450=6.4 mJy. 


30 min scans reduced using the dynamic iterative map-maker of 
the SMURF package (Jenness et al. 2011), flat-fielding, and then 
solving for model components assumed to make up the bolome¬ 
ter signals (atmospheric, astronomical, and noise terms). The sig- 


2.2 IRAM/PdBI Observations 

The IRAM-PdBI observations were taken in D configuration us¬ 
ing the 5-antenna sub-array (properties of the observing runs and 
resulting synthesized beam sizes are summarized in Table 1), at a 
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pointing center: RA=17 01 17.779, Dec=+64 14 37.85, J2000.0, 
taken from the SCUBA-2 centroid. Flux calibration was achieved 
by observing various calibrators (3C273, 3C345, B0234+285, 
B1749+096). The quasar B1300+580 was used as a phase and am¬ 
plitude calibrator. Data were processed using the most recent ver¬ 
sion of the GILDAS software. We resampled the cubes in 90 km 
s” 1 channels, and imaged them using the GILDAS suite mapping, 
adopting natural weighting. The full PdBI spectrum is shown in 
Fig. 2. To obtain flux measurements, we deconvolved the visibili¬ 
ties using the CLEAN task with natural weighting, and applied the 
corresponding primary beam correction (negligible at the small 3” 
offset of the F1S1700.850.1 position). 

2.3 SMA observations 

Followup SMA observations were performed on August 15 and 
September 7. 2014 in the compact and extended configurations 
respectively (beam sizes ~2" and ~0.85” respectively) in good 
weather (T 225 GHZ ~ 0.08) with a total on-source integration time 
of approximately 6 hr. The USB was tuned to 345 GFlz, and com¬ 
bined with the LSB for an effective bandwidth of ~ 4 GHz at 340 
GHz, which yielded a final rms of 1.1 mjy. The pointing centre 
was the same as that for the PdBI observations described above. 
The data were calibrated using the MIR software package (Scov- 
ille 1993), modified for the SMA. Passband calibration was done 
using 3C 84, 3C 111, and Callisto. The absolute flux scale was set 
using observations of Callisto and is estimated to be accurate to 
better than 20%. Time-dependent complex gain calibration was 
done using 1858+655 (0.6 Jy, 21.8° away) and 1827+390 (1.8 Jy, 
37.7° away). We detect HS1700.850.1 in the synthesized image at 
~ 13a. The calibrated visibilities were best fit by a single point 
source with an integrated flux density of Sgg o/im = 14.5 ± 1.1 mjy 
at a position of a(J2000) =17 01 17.779 and <)(J2000) = +64 14 
37.85, consistent within 0.3” with the centroid of the line and con¬ 
tinuum measurements from PdBI. The astrometric uncertainties are 
Aa = 0.24arcsec (0.20arcsec systematic; 0.13arcsec statistical) 
and AS — 0.22arcsec (0.19arcsec systematic; 0. lOarcsec statisti¬ 
cal). 

3 RESULTS 

3.1 The redshift search 

Our original IRAM-PdBI observations searched for the 12 CO and 
[Cl] lines in HS 1700.850.1 at a redshift ~ 2.30. This was moti¬ 
vated by it lying near the core of a massive z = 2.31 protocluster 
in the Keck Baryonic Structure Survey (Steidel et al. 2011), and 
having a plausible 4.5^m identification with a complex of known 
z = 2.305 BX galaxies. Being such an extreme source we would 
not expect to find it by chance in the small 0.03 deg 2 H1700+64 
pointing, or even in the 0.19 deg 2 survey of the six KBSS fields 
without the high density environment of the proto-cluster in this 
field. Using typical single dish 850/rm counts (Coppin et al. 2006; 
Weiss et al. 2009) we expect a surface density of just 0.01 deg” 2 
for such 19 mjy sources. The source counts of the brightest submm 
sources (S870^m>9 mJy) measured by ALMA fall more steeply 
than those measured with lower resolution single-dish observations 
due to the multiplicity of the brightest sources when observed at 
arcsecond resolutions (e.g., Karim et al. 2013). As HS1700.850.1 
retains at least ~15mJy in the interferometric followup (and as 
much as 19mJy if the compact configuration SMA data is taken 


at face value), the expected >15mJy source density drops by an¬ 
other factor of ten based on the ALMA-derived count of Karim et 
al. (2013) or Simpson et al. (2015b). 

The IRAM-PdBI observations showed no such lines at the ex¬ 
pected frequencies, and thereby implied that the proposed identifi¬ 
cation and redshift of z ~ 2.31 was incorrect. As the 7.2 GHz of 
frequency coverage did reveal a strong line detected serendipitously 
at 150.6 GHz near the edge of the 2mm setting, we began a search 
to identify the correct redshift. The centroid of this 150.6 GHz line 
detection, as well as the PdBI 2mm and 3mm continuum detections 
of HS 1700.850.1, both localize the emission from the source within 
optical/IR imagery to within ~ 0.5” (Fig. 1), while followup SMA 
870/rm continuum imaging unequivocally associates the 850/rm 
SCUBA-2 source to a single IRAC galaxy which is well detected in 
the /7STF814W image (Fig. 1) revealing a faint disturbed, tadpole¬ 
like galaxy. Ground-based imagery marginally detects the galaxy at 
7Z = 26.1, g = 27, U > 27.1, making optical spectroscopy very 
difficult, and the z = 2.82 redshift would be difficult for near-IR 
spectroscopy, although the R' 3 ,ab = 23.2 is feasible in general 
with a 10m telescope. Note that these colors (e.g., g' — 1Z = 0.9 
mag) make sense for a heavily reddened object at z= 2.82. However, 
the detection at (/-band, places loose constraints on the redshift to 
z < 4, while the comparable IRAC fluxes, S(4.5/.im)=10.2/Hy .and 
S(8.0/im)=13.2/iJy suggests the 1.6pm stellar bump should lie be¬ 
tween, ranging between 2.4<z<3.7. The very red K s — 4.5pm 
coloyr (1.8 mags AB) is large, though consistent with a high stel¬ 
lar mass object. The MIPS 24/jm flux of 0.17mJy compared to 
the S 85 o M m flux suggests that PAH emission is not contributing 
to the 24pm band, pushing the redshift to z > 2.5, as illustrated 
in the SED template fit (Fig. 4). Our final constraint on the red¬ 
shift comes from the long-wavelength photometry (Fig. 4, Table 2), 
where 3mm, 2mm, 870/rm, and 450^m fluxes suggest a similar 
range of 2.5<z<4.5 for a likely range of SMG dust temperatures, 
25 K<T d <45 K. 

While these constraints may not seem very useful, the fact 
that we have a solid line detection at 150.6 GHz means that only 
three possible 12 CO identifications provide redshifts in this range: 
12 CO(4-3) at z = 2.82, 12 CO(5-4) at z = 3.58, or 12 CO(6-5) at 
z = 4.37. Aided by these photometric redshift estimates, in two 
subsequent PdBI followup observations we searched for z = 3.58 
12 CO(4-3) at 100.4 GHz, and then at 91.5 GHz searched simul¬ 
taneously for z = 2.82 12 CO(4-3) and [Cl] at z = 4.37. With 
a detection of a second line (Fig. 2) we confirmed the redshift 
z=2.816 for one of the most intrinsically luminous galaxies in the 
distant Universe (Lfir ~ 5 x 10 13 Lq, § 3.3) and secured a sec¬ 
ond 12 CO line for astrophysical diagnostics. Blind 12 CO redshifts 
for unlensed SMGs are still exceedingly rare in the literature. We 
are only aware of Daddi et al. (2009) GN20 at z = 4.05 (a truly 
serendipitous discovery), and the z = 5.2 for HDF850.1 (Walter et 
al. 2012). 

3.2 Source properties 

HS1700.850.1 is very bright at 850/rm. but shows no sign of 
lensing - the nearest galaxy detected is the low-mass star form¬ 
ing galaxy BX951 2.9" away (z=2.3053, M* = 2 x 10 9 M 0 ). 
HS 1700.850.1 is thus not suffering from differential lensing and as¬ 
sociated modelling issues and biases that could affect our interpre¬ 
tation of the 12 CO lines. The 12 CO spectra are shown zoomed in at 
a velocity resolution of 60 km s” 1 in Fig. 3. Both lines are detected 
at high significance (11.3 a and 6.9 a for continuum subtracted (3- 
2) and (5-4) lines respectively, in integrated intensity). The line 
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Figure 4. The SED of HS1700.850.1, with the Swinbank et al. (2014) com¬ 
posite SMG template overlaid. Significantly hotter or cooler SEDs than this 
~ 35 K template are reasonably rejected by the photometry. The template 
is split at a rest frame lOfim in order to fit the optical/near-IR faint compo¬ 
nent as decoupled from the luminous far-IR to submm dust component. The 
far-IR luminosity of this model is 3.8 • 10 13 Lq. 
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58.53 GHz for the 12 CO line-pairs at redshifts 2.82 and 6.64 re¬ 
spectively. Our observations yield 6 u = 58.58 ± 0.02 GHz, clearly 
consistent with the former, and supported by the photometric red- 
shift. 

With the precise redshift and the observed 12 CO line lumi¬ 
nosities in hand, we can estimate the molecular gas content of 
HS1700.850.1. The observed global 12 CO(5-4) to 12 CO(3-2) line 
ratio (0.67±0.09) implies that the CO emission is sub-thermally 
excited, at least for the CO(5^f) line. The average line ratio re¬ 
ported by Bothwell et al. (2013) for a sub-sample of their 30 
SMGs is r53=0.62. This line ratio is similar to that observed for 
SMM J16359+6612 (Weiss et al. 2005) and we use a similar scal¬ 
ing to estimate a 12 CO(1-0) line luminosity of L' « 1.86 • 
10 11 K kms -1 pc 2 . This translates into a molecular gas mass of 
Mh 2 ~ 2.1 ■ 10 11 Mq using a standard ULIRG conversion fac¬ 
tor of 1.0M & (K kms' 1 pc 2 ) -1 found for typical 2 ~ 2.5 SMGs 
in Bothwell et al. (2013), and including a 15% contribution from 
Helium. (Downes et al. 1998 adopted 0.8 Mq (K kms -1 pc 2 ) -1 ). 


Table 3. Gaussian decomposition of the CO line spectra. Values for the total 
line are estimated from integrating the continuum-subtracted spectra, while 
the red and blue component properties are estimated from the Gaussian fits. 
Errors in redshift estimation are d2=0.001. The measured FWHM for the 
total line is quoted as a, 2.355 X smaller. L'CO is in units of 10 10 K km 
s -1 pc 2 . Note that for ct(CO)=l Mq (K km s -1 pc 2 ) -1 , as estimated in 
Bothwell et al. (2013), L'COi—o - M(gas), which should then be corrected 
by an additional 15% for Helium. 



blue 

red 

total 

2 (3—2) 

2.812 

2.820 

2.816 

5v(3-2) [kms” 1 ]" 

-305 ± 9 

±280 ± 11 

0 

5v(5-4) [kms” 1 ] 0 

-280 ±21 

±250 ± 19 

0 

cr(3-2) [kms -1 ] 

140 ± 25 

142 ± 20 

400 ± 57 c 

cr(5-4) [kms -1 ] 

200 ± 45 

143 ± 24 

349 ± 45 c 

I 3-2 [Jy km s -1 ] 

1.26 ±0.16 

0.88 ±0.13 

2.32 ±0.21 

h -4 [Jy km s” 1 ] 

1.84 ± 0.23 

1.5 ± 0.21 

3.44 ± 0.30 

L'COs -2 

5.18 

3.62 

9.65 

L'C0 5 -a 

2.72 

2.21 

5.08 

L'COi-o b 

9.95 

6.95 

18.56 


a Center velocity relative to z— 2.816. 

b Derived using average i' 3 i = 0.32 scaling relation of Bothwell et al. 
(2013). 

c These are the best fit single Gaussian to the entire line profile. 

profiles for both lines are similar and well described by a double 
Gaussian with a peak separations of 585 and 530 kms -1 for (3-2) 
and (5^1) 12 CO lines respectively, and individual peak FWHMs of 
~330-470 kms -1 . While the (3-2) line shows an apparent peak 
flux asymmetry compared with the (5-4) peaks, in fact the flux ra¬ 
tios of the blue and red lines are quite similar, r53^;„ e =0.69±0.12 
and r53 j7 . e( i=0.59±0.09. The parameters derived from Gaussian fits 
to both line profiles are given in Table 3. We find only marginal 
evidence for any spatial separation in the two peaks of either of the 
12 CO lines < 1". The frequencies unambiguously identify the lines 
as CO(3-2) and CO(5^1). Combining the centroids of both lines, 
we derive a variance-weighted mean redshift for HS1700.850.1 of 
2 = 2.816±0.001. While the observed frequencies might also be 
interpreted as 12 CO(6-5) and 12 CO(10-9) at 2 = 6.64, the 12 CO 
ladder is not equidistant in frequency which results in small differ¬ 
ences for the frequency separation of the line-pairs as a function of 
rotational quantum number. The frequency separation is 58.58 and 


3.3 SED modelling 

The SED of HS1700.850.1, with the Swinbank et al. (2014) com¬ 
posite SMG template overlaid is shown in Fig. 4. Significantly hot¬ 
ter or cooler SEDs than this ~ 35 K template are reasonably re¬ 
jected by the photometry. The template is split at a rest frame 10/rm 
in order to fit the optical/near-IR faint component as decoupled 
from the luminous far-IR to submm dust component. The far-IR 
luminosity of this model is 3.8 x 10 13 Lq, translating to a SFR of 
2500 M 0 yr -1 (Kennicutt 1983). HS1700.850.1 is not unusually 
faint in optical/IR magnitudes for SMGs with Ss5o M m~3-5 mJy, 
but its far-IR output is significantly higher, and it is thus propor¬ 
tionally far more obscured than ALESS SMGs (see Swinbank et al. 
2014). 

A hyper -2 fit to the IRAC and near-IR photometry yields a 
stellar mass of 3X10 11 Mg, assuming a Chabrier IMF. Together 
with the above gas mass estimate, the gas fraction is ~40%, and 
precludes using substantially higher conversion factors, a(CO)^> 
1 , which would drive the gas fraction up towards 100%. 

3.4 SMA and PdBI continuum maps and source multiplicity 

The PdBI observations reveal only a single significant detection 
in both continuum and line emission, showing that there is no 
obvious multiplicity at 2mm/3mm or in CO lines. As such there 
are no nearby (< 25" primary beam) 12 CO-luminous compan¬ 
ions to HS 1700.850.1. Similarly the SMA image does not de¬ 
tect any additional significant sources. The combined 3mm con¬ 
tinuum map reaches 0.35 mJy rms. At this depth we are sensitive 
to SMGs with SFRs~1000Mq yr -1 . The 2mm continuum map 
reaches 1.2 mJy rms and should detect similar SFRs for a typical 
/3 ~ 1.5 dust SED. The SMA 870/rm map rms of 1.1 mJy is much 
more sensitive to dusty sources, and should detect galaxies with 
SFRs~3OOM 0 yr -1 . We note that a CO-luminous DRG in this 
field with a SFRs~200 M 0 yr -1 from Chapman et al. (2015) is not 
detected in any of these continuum maps, nor are any of the other 
known galaxies lying within the primary beam. 

By considering only the SMA compact configuration data (2" 
beam size), we measure a larger flux density of 19.3±2.8 mJy, close 
to that measured by SCUBA-2 in the 14.5" beam, although the 
large measurement error agrees with the extended configuration 
flux measurement within ~2<r. Nonetheless, this may suggest there 
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is either a modest conflict in the flux scale between compact and 
extended configurations, or that the source may be marginally re¬ 
solved on scales greater than the average beam size achieved in our 
extended configuration data, > 0.85”. However, this is unusual for 
bright SMGs (Simpson et al. 2015a) which typically are resolved 
with ~ 0.3 — 0.4 source sizes. Ubiquitous SMG multiplicity results 
from surveys (e.g., Wang et al. 2010, Hodge et al. 2013, Simpson 
et al. 2015b) and theoretical arguments (e.g. Hayward et al. 2012, 
2013a, 2013b) would argue that there may well be faint compan¬ 
ions within the SCUBA-2 beam, and the similarity of the compact 
SMA data flux may just be a lucky coincidence. 


4 DISCUSSION 
4.1 Gas dynamics 

At first glance the double peaked line profile suggests a circum- 
nuclear molecular toroid. However, if such a stable molecular 
gas distribution would supply the large star formation rate (« 
2500 M 0 yr _1 ), we would expect a high gas excitation. The ob¬ 
served line ratio rs3=0.67±0.09 is low compared to well-studied 
nearby starburst galaxies NGC 253 and M 82, although it is typical 
of that measured by Bothwell et al. (2013), a brightness temper¬ 
ature ratio r53=0.62, for SMGs with characteristic Lssop?™ lower 
than HS 1700.850.1 by a factors of 3x on average. However, the 
two line peaks could also be two disk components merging, indi¬ 
vidually unresolved. The line integrated gas properties between the 
spectral components are similar in the red and blue component line 
ratios (r53, re d=0.6±0.1, r53,biue=0.7±0.2), contrasting the compo¬ 
nent line widths and central velocities which do differ between the 
different excitation lines. This argues that the (5-4) line may be 
coming from a different volume of gas than the (3-2) line, which 
could be true in either a disk-like or merger configuration. If the 
SMA dust continuum size is the same as the CO gas (< 0.85 " or 
< 6.7 kpc), then either the disk or the merger is in a very com¬ 
pact state. We note that a ring-like rotating disk distribution of gas 
was strongly endorsed as an explanation for another similar double- 
peaked SMG, SMM J02399-0136 (Genzel et al. 2003), however 
from multi-frequency, high-resolution line mapping observations, 
Ivison et al. (2010) showed that the system comprises a multiple 
merger. 

HS1700.850.1 is clearly a very massive system, and is in 
particular a very massive barvonic system. Observational stud¬ 
ies of local ellipticals (Keeton 2001; Boriello, Salucci, & Danese 
2003; cf. Loewenstein & Mushotzky 2003) suggest that typically 
~ 30% of the total mass within the optical effective radius will be 
dark; the interior regions of massive local disk galaxies have sim¬ 
ilarly low dark matter fractions (Combes 2002). By analogy with 
these local systems, HS 1700.850.1 should have a baryonic mass 
of ^ 4 x 10 11 sin -2 i~ x M 0 within 7 kpc; M gaa ~ 2 x 10 11 Mg 
then implies that a limit to the stellar mass within the same radius is 
~ 2 x 10 11 Mg. The Kennicutt (1983) stellar initial mass function 
(IMF) models of Cole et al. (2001) imply that M* ~7x 10 10 M 0 
for the galaxy stellar mass function in the local Universe. Based on 
the stellar content of its inner 7 kpc, then, HS 1700.850.1 is already 
at least a 3M* system, on its way to becoming a ^ 4M* system as 
it turns the remainder of its molecular gas into stars. 



L'C0(3-2) K km/s pc" 2 

Figure 5. SFR versus L'CO(3-2) for various literature galaxies compared 
with HS1700.850.1 (normal star forming galaxies - Tacconi et al. 2013; 
SMGs - Bothwell et al. 2013). CO(3-2) is adopted for maintaining all 
displayed galaxies as close to observed transitions as possible. We show 
L'(CO) for HS1700.850.1, revealing it is nominally under-luminous in 
12 CO compared to other ultra-luminous normal star forming galaxies, but 
not unusual compared to many SMGs in Bothwell et al. (2013) in lying a 
factor 2-3 X above the Tacconi et al. (2013) best fit relation. 


4.2 An Eddington-limited starburst? 

If we assume that HS 1700.850.1 is starburst dominated and make 
assumptions about the morphology and kinematics, we find it may 
be radiating close to or at the Eddington limit of its starburst. 
Elmegreen et al. (1999) demonstrated that feedback from ongo¬ 
ing star formation sets a physical limit on the minimum size of 
a star forming region. Radiation pressure from high luminosity star 
formation regions should produce strong momentum-driven winds 
due to high dust opacity seen by the ultraviolet light produced by 
young stars (e.g., Netzer & Elitzur 1993). These winds are confined 
by the gravitational potential, which scales as <f> ~ f g a 2 log D for 
an isothermal sphere (with the gas fraction f g , stellar velocity dis¬ 
persion a, diameter of the starburst region D). In the optically thin 
limit (appropriate for optically thick clouds with a small volume 
filling factor embedded in the diffuse interstellar medium) and as¬ 
suming a Chabrier et al. (2003) initial mass function (IMF) this 
leads to a maximum star formation rate 

SFR max = 360 ufooOkpcKioo M 0 yr -1 (1) 

We measure the scale size of the starburst _Dk pc as the Gaussian 
FWHM of the SMA detected source, and <7400 as the line-of-sight 
12 CO gas velocity dispersion in units of 400 km s” 1 which is the 
best fit single Gaussian to the entire CO(3-2) line ; k 100 is the 
dust opacity in units of 100 cm 2 g _1 , which we set to 1, but note 
that many dust models allow for significantly higher opacity, and 
in particular for the ultraviolet radiation produced by young mas¬ 
sive stars during a starburst (e.g., Li et al. 2001). HS 1700.850.1 
has a very high far-IR luminosity, with Lfir « 4 x 10 13 L 0 
and SFR « 2500 M 0 yr -1 on a physical scale of < 7 kpc. This 
very luminous SMG is therefore near the Eddington limit ^ 2000 
M 0 yr -1 for a starburst on those scales. The size/SFR ratio of 
HS1700.850.1 is typical of the more luminous SMG components 
identified with ALMA observations by Simpson et al. (2015a), sug¬ 
gesting near Eddington-limited starbursts are not unusual in the 
SMG population. This motivates obtaining additional higher res¬ 
olution 870/rm and 12 CO observations of HS1700.850.1 to truly 
resolve the emission of this bright galaxy. 
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4.3 Lqo to Lfir comparison with literature 

We next turn to a discussion of HS1700.850.1 versus other 12 CO- 
detected sources, plotting Leo -SFR for the known and published 
high-z sources (Fig. 5), and translating Lfir to SFR in the SMG 
populations using standard conversions (Kennicutt 1983). Of note 
is that FiS 1700.850.1 has a significantly higher Lgo than any of the 
35 unlensed SMGs in the Bothwell et al. (2013) survey. Compara¬ 
bly 12 CO-luminous SMGs were found in a few cases - Ivison et al. 
(2013) identified a binary HyLIRG merger from a ~200 deg 2 Her- 
schel survey, while one lensed SMG from the ~1300 deg 2 South 
Pole Telescope survey appears to be this luminous in 12 CO after 
a demagnification factor of ~ 5x (Hezaveh et al. 2013). Clearly 
finding SMGs with such high h' co has not been an easy task, even 
probing the widest field mm-wave surveys. 

F1S1700.850.1 is somewhat under-luminous in 12 CO for its 
SFR (or a high efficiency star former) compared to the average re¬ 
lation from Tacconi et al. (2013) for normal star forming galaxies 
(shown in Fig. 5). Flowever, compared to SMGs front Bothwell et 
al. (2013), F1S1700.850.1 is not unusual in lying 2-3x above the re¬ 
lation, an observation that supports application of the ULIRG con¬ 
version rate from Leo to total molecular gas mass. 


4.4 The environment of HS1700.850.1 

Given the extreme nature of HS 1700.850.1, we explore whether 
the surroundings of this SMG are unusually clustered. The envi¬ 
ronment of HS 1700.850.1 is traced in great detail by the LBG sur¬ 
vey in the field, presented in Fig. 6 as a histogram binned to 500 km 
s" 1 to highlight any possible neighbouring galaxies. HS1700.850.1 
appears to reside in a sparse region, with the closest (and only) 
neighbouring galaxy within ±1500 km s^ 1 at z=2.807 and lying 
5' (2.5 Mpc) away. The LBG selection function appropriate to this 
field (see Steidel et al. 2014) suggests that only two LBGs are ex¬ 
pected in a 1500 km s _1 region from the 243 spectroscopic LBGs 
with measured redshifts in this field, thus finding only one about 
HS1700.850.1 is not statistically significant. There is one modest 
over density of six sources at z = 2.751, which includes the QSO 
HS1700+64 itself, contrasting with the sparse HS1700.850.1 envi¬ 
ronment nearby. Off the figure, at z = 2.31, there is a large over¬ 
density consistent with a forming proto-cluster (Steidel et al. 2005). 

Using the extensive spectroscopic database of the KBSS sur¬ 
vey, Trainor & Steidel (2012) study similar extreme sources to 
HS1700.850.1 at z=2.5-3.5 - HyLIRG quasars (HLQSOs). They 
find the HLQSOs themselves are associated with a 7 overden¬ 
sity in redshift when considered on scales of ~5 ft -1 Mpc, corre¬ 
sponding to a velocity scale of ov, P ec = 200 km s^ 1 and a pro¬ 
jected scale of R~200 proper kpc. By studying the redshifts of 
galaxies around each HLQSO, they find however that there are no 
overdense peaks of similar significance. Repeating this procedure 
on random galaxy redshifts shows that the HLQSOs are correlated 
with much more significant small-scale overdensities than the aver¬ 
age galaxy. They deduce that each HLQSO in their sample inhabits 
a DM halo with mass log(M h /Mg) > 12.1±0.5, or a median halo 
mass of logtMh.med/Mg) = 12.3±0.5. Crucially, the number den¬ 
sity of these halos exceeds the number density of HLQSOs by a 
factor ~ 10 6 — 10', prompting the conclusion that the host halos 
of HLQSOs are not rare. Clearly HS1700.850.1 follows this trend, 
inhabiting an environment that is indistinguishable from the aver¬ 
age field galaxy density in Fig. 6. 

The scarcity of HLQSOs is likely due to an extremely improb¬ 
able small-scale phenomenon that produces HLQSOs, possibly re- 



2.6 2.8 3 

z 

Figure 6. The environment of HS1700.850.1 derived from the LBG survey 
in the field, with 500 km s —1 bins to highlight any possible neighbouring 
galaxies (the closest one at z=2.807 lies 5' away). HS 1700.850.1 appears to 
reside in a relative void, although the selection function suggests only two 
expected LBGs in a 1500 km s —1 region for the 243 spectroscopic galax¬ 
ies from z = 1.5-3.5 measured in this field (the overall redshift selection 
function for the continuum-selected redshift survey in all 15 fields of KBSS 
is shown as a blue line - see e.g. Steidel et al 2014). There is a 6-object 
galaxy spike at the redshift of the QSO HS1700+64 itself, contrasting the 
sparse environment of HS1700.850.1, while at z = 2.31 (off the plot axes) 
there is in fact a very large galaxy over-density consistent with a forming 
proto-cluster (Steidel et al. 2005). 


lated to an atypical galaxy interaction geometry or similar scenario, 
or else extremely short durations of the events. Analogously in 
HS1700.850.1 and other extreme SMGs, an unlikely merger trajec¬ 
tory perhaps with unusually gas-rich progenitors (e.g., Narayanan 
et al. 2010) may explain rare bursts approaching Ssso/rm ~20 mJy 
and Lir ~ 2 x 10 13 Lq (e.g., Ivison et al. 2013). This mo¬ 
tivates obtaining more detailed followup of the resolved double- 
peaked 12 CO emission in HS 1700.850.1 to better understand the 
circumstances generating such enormous star formation rates and 
illuminated dust masses. We might also hypothesize that a very 
short duration starburst could explain the extreme luminosity of 
HS1700.850.1. For example, if we assumed that every Sg5o Mm =15- 
20 mJy SMG was a brief episode in the more numerous 5 mJy SMG 
population (1000X more numerous in the ALMA-based counts of 
Simpson et al. 2015b, but still very luminous starbursts in them¬ 
selves), then the HyLIRG burst would last only 50kyr of a charac¬ 
teristic burst of 50Myr (e.g., Chapman et al. 2005, Tacconi et al. 
2006). This may not be a well motivated scenario - such a short 
timescale could likely only be driven by intense, hot AGN activ¬ 
ity, yet there is no sign of an AGN in HS1700.850.1. Further, an 
AGN could not boost the 850/rm emission while keeping the typi¬ 
cal SMG appearance in Fig. 4. Most likely the SMG HyLIRGs are 
driven by a combination of factors. 

Another estimate on the halo mass can be obtained using 
dynamical mass estimates from the peculiar velocities and the 
projected scale of the galaxy overdensities. With this approach 
Trainor & Steidel (2012) estimate HLQSOs are on average asso¬ 
ciated with even larger group-sized environments with total mass 
log(M grp /M©)~13. The SMG clustering analysis in Hickox et al. 
(2012) suggested that the SMG population is hosted by haloes of 
comparable mass (see also Blain et al. 2004). This overdense envi¬ 
ronment with small relative velocities would increase the probabil¬ 
ity of a rare event, but an unusual merger configuration is still likely 
required to generate such large luminosities. In HS1700.850.1, 
even a group environment of this mass is unlikely based on the 
lack of explicit companions within 1500 km s _1 and within ~2.5 
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Mpc. However, a ~ 10 13 Mg halo will typically only have a cou¬ 
ple of bright LBG-class galaxies in it, and might be easy to miss in 
an incomplete spectroscopic survey. 

The Trainor & Steidel (2012) analysis has clearly demon¬ 
strated that HyLIRG-class sources do not require environments 
very different from their much less luminous counterparts. Ev¬ 
idently, the exceedingly low space density of at least HLQSOs 
(~10 -9 Mpc -3 ) results from a one-in-a-million event on scales 
<C 1 Mpc, and not from being hosted by rare dark matter ha¬ 
los. HS 1700.850.1 appears to inhabit a similarly typical halo. A 
moderate mass environment was also noted for a highly signifi¬ 
cant association of 8 SMGs at z~1.99 by Chapman et al. (2009), 
whilst revealing that z=2-3.5 SMGs generally do not align well 
with LBG overdensities within the GOODS-N survey field. Similar 
anti-biasing of SMGs is noted in the ECDFS field (A. Danielson in 
preparation). Notable exceptions exist at higher redshifts, including 
a z=4.05 galaxy overdensity found through an association of bright 
SMGs (Daddi et al. 2009), the COSMOS z=5.1 structure (Capak et 
al. 2011), and HDF850.1, which apparently lies in a z=5.3 galaxy 
overdensity (Walter et al. 2012). 

Miller et al. (2015) used the Hayward et al. (2013a) model for 
the SMG population to estimate typical environments of the bright¬ 
est SMGs and associations of multiple SMGs found in ~500 Mpc 3 
volumes. This model is based on mock galaxy catalogues derived 
from the Bolshoi simulation (Klypin et al. 2011; Behroozi et al. 
2013) and insights obtained from performing radiative transfer on 
hydrodynamical simulations (Hayward et al. 2011; Hayward et al. 
2013b). Miller et al. (2015) find the SMGs inhabit a range of mass 
scales going down to the average cosmic density in the simulation, 
and rarely are found in the most massive regions at these z = 2- 
3 epochs (although at significantly higher redshifts they are better 
tracers of the most massive regions). The surrounding LBG sur¬ 
vey indeed suggests HS1700.850.1 is statistically consistent with 
inhabiting the average cosmic density. 


5 CONCLUSIONS 

We presented an IRAM-PdBI spectral survey of the HyLIRG 
HS1700.850.1, the brightest 850/.tm source found in the SCUBA-2 
survey of the HS 1700+64 field, and indeed all six SCUBA-2/KBSS 
fields. Our primary goal was to identify its redshift, guided hy pho¬ 
tometric redshift estimates in the optical through far-IR, and study 
its properties and environment. We conclude: 

• The redshift for HS1700.850.1 is z=2.816, secured through 
the 12 CO(3-2) and (5-4) lines, validating the use of crude photo¬ 
metric redshifts to guide 12 CO redshift searches. 

• The clear doubled-peaked profile of HS 1700.850.1 in both 
CO transitions is likely a result of an end stage merger, within 6 kpc, 
as constrained by the size of the 850/im emission. 

• HS 1700.850.1 exhibits typical Lq D for its far-IR luminosity 
or SFR compared to other unlensed SMGs, although its Lq 0 is 
significantly larger than any SMGs identified in the Bothwell et al. 
(2013) survey. 

• The environment of HS 1700.850.1 is indistinguishable from 
the average field galaxy density, and no nearby LBG companions 
are identified within 1500 km s -1 . The HyLIRG burst could result 
from a 10 -6 unlikely event on scales <§C 1 Mpc, and not from being 
hosted by a rare dark matter halo, but could also be due to a very 
short ~ 50 kyr burst timescale. 
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